Structural signatures of mobility on intermediate time scales in a supercooled fluid 
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We use computer simulations to explore the manner in which the particle displacements on in- 
termediate time scales in supercooled fluids correlate to their dynamic structural environment. The 
fluid we study, a binary mixture of hard spheres, exhibits classic signatures of dynamic heterogene- 
ity, including a bifurcated single-particle displacement distribution (i.e., subpopulations of immobile 
and mobile particles). We find that immobile particles, during the course of their displacements, 
exhibit stronger average pair correlations to their neighbors than mobile particles, but not neces- 
sarily higher average coordination numbers. We discuss how the correlation between structure and 
single-particle dynamics depends on observation time. 



When fluids are supercooled (or overcompressed) to- 
ward their glass transition, their single-particle dynam- 
ics undergo qualitative changes*^— One example is the 
emergence of dynamic heterogeneity at time scales inter- 
mediate between ballistic and diffusive regimes of par- 
ticle motion, which manifests as spatially dependendent 
relaxation processes in the liquid 2 -! 6 -^— and bifurcated 
(multi-peaked) probability distributions associated with 
single-particle displacements. The latter suggests the 
presence of distinct subpopulations of particles with dif- 
ferent mobilities on these time scalesii£~— Dynamic het- 
erogeneities of this sort continue to attract wide interest 
because of their perceived consequences for other pro- 
cesses in deeply supercooled liquids including the break- 
down of the Stokes-Einstein relationship 6 ' 18 ^ 3 " — and the 
emergence of non-exponential trends in relaxation of the 
structure factori 6 ^ 

Several recent studies have systematically explored the 
extent to which a fluid's particle configuration at a given 
time influences the spatial distribution of relaxation pro- 
cesses that immediately follow i^Zr— Other investigations 
have focused on probing how certain static structural 
properties of fluids can be used to correlate the effects 
that temperature, density, external fields, and interpar- 
ticle potential have on long-time dynamics.— ~— Here, we 
explore a different, but related, question. Does the local 
structure surrounding particles in a deeply supercooled 
liquid, averaged over an intermediate time scale relevant 
for dynamic heterogeneities, strongly correlate with their 
mean-square displacements during the period of observa- 
tion? In other words, do particles with high mobility 
sample more disordered structural environments during 
the course of their displacements than those with low 
mobility and vice versa? 

We investigate this question via event- driven^ molec- 
ular dynamics simulation of dense binary fluid mixtures 
of hard spheres. We set the ratio of particle diam- 
eters in these mixtures to a\j<y% = 1.3 and the ra- 
tio of particle masses to mxjm^ = (ci/^) 3 , parame- 
ters that mimic concentrated colloidal suspensions that 
were recently investigated experimentally— We simulate 
Ni = N 2 = 1000 particles in a periodically-replicated 
cubic cell of volume V. We present results for parti- 



cle packing fractions <p = ir(Niaf + N 2 <J2)/6V of 0.57 
and 0.582, which, as we show below, correspond to state 
points with unimodal and bifurcated displacement dis- 
tributions on intermediate time scales, respectively. For 
brevity, we report quantities that are implicitly nondi- 
mensionalized by appropriate combin ations of the length 
scale, l c = a-2 and time scale t c — y ^erl/fceT, wnere 
fee is Boltzmann's constant. We focus on the dynam- 
ics and structure of the smaller type 2 particles, but we 
note that the behavior of the larger type 1 particles (not 
shown here) is qualitatively similar, as might be expected 
given the mild particle-size asymmetry of the fluid. 




FIG. 1: (Color online) Mean-square displacement (Sr 2 ) of the 
smaller (type 2) particles versus time t at packing fraction 
ip — 0.57 (black) and 0.582 (red). Dashed lines are fits of 
the Einstein relation (Sr 2 ) = 6Dt to the long-time behavior, 
resulting in tracer diffusion coefficients of D = 1.2 x 10 -3 and 
1.5 x 10~ 4 at (p = 0.57 and 0.582, respectively. 

We begin by examining the time t dependence of the 
average mean-square displacement, (Sr 2 ), for the type 
2 particles. In particular, Figure [T] displays results for 
packing fractions of ip = 0.57 and 0.582. At both state 
points, the fluid exhibits a mean-square displacement 
plateau at intermediate times, which is characteristic of 
"cage" dynamics^ 3 Schematically, the plateau separates 
the ballistic motion that occurs at very short times (be- 
fore motion is temporarily hindered by collisions with 
the cage of nearest-neighbor particles) and the diffusive 
motion that particles ultimately exhibit at long times (af- 
ter breaking through the cage). As should be expected, 
the plateau occurs at smaller displacements and persists 
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FIG. 2: (Color online) Distribution of the logarithm of mean- 
square displacement P[log 10 (<5r 2 )] for type 2 particles of the 
binary hard-sphere mixture described in the text. Data cover 
several time intervals t — r for packing fractions (a) if — 
0.57 and (b) <p = 0.582. The numbers in the figure provide 
the value of r for each curve. In (b), the vertical line is at 
Sr 2 — 0.56, the approximate dividing line between mobile and 
immobile particles discussed in the text. 



for longer times when tp is increased, indicating that the 
cage formed by the nearest-neighbor coordination shell 
becomes both tighter and more difficult to disrupt at 
higher packing fraction. 

In order to characterize dynamic heterogeneities of the 
type 2 particles at the state points examined above, we 
follow Ref . and investigate the distributions of the log- 
arithm of the mean-square displacement (P[log 10 (5r 2 )]) 
over a variety of time intervals (t = r) . Specifically, Fig- 
ure [2a| displays P[log 10 (<5r 2 )] for the fluid at <p = 0.57 
and values of r that correspond to average mean-square 
displacements which span from the plateau "cage" re- 
gion to the beginning of the diffusive regime (see Fig- 
ure [T]). The main point of Figure [5a] is that, for all 
values of t, P[\og 10 (Sr 2 )] remains unimodal, indicating 
that pronounced single-particle dynamic heterogeneities 
have not yet emerged in the fluid at this packing fraction. 
However, also notice that for the displacement distribu- 



tions corresponding to intermediate times, particularly at 
t = 50, a shoulder at higher mean-square displacements 
becomes evident. This shoulder is a precursor to the bi- 
furcated displacement behavior that occurs at cp = 0.582, 
which we discuss in detail below. 

Figure [2b] displays the behavior of P [log 1Q (Sr 2 )] at 
ip = 0.582. Notice that, for intermediate times be- 
fore diffusive behavior is reached (500 < r < 4000), 
P[log 10 (Sr 2 )] shows signatures of bimodality, i.e., two- 
peak structure. This suggests that subpopulations of 
mobile and immobile particles are emerging. A quali- 
tative dividing line between the two subpopulations can 
be drawn at 8r 2 « 0.56 (vertical dashed line in FigurcQ]), 
which distinguishes the particles that are still caged (i.e., 
localized) on intermediate time scales from those that 
have broken through their nearest-neighbor coordination 
shells to attain larger displacements. For descriptive pur- 
poses, we label particles with Sr 2 < 0.56 "immobile" (on 
these time scales), while we label those with Sr 2 > 0.56 
"mobile" . The type of dynamic behavior depicted in Fig- 
ure [2b] has been well documented in other systems^ Be- 
low we investigate whether the immobile particles experi- 
ence, on average, a more ordered structural environment 
than the mobile particles during the time intervals of 
their respective displacements. 

To carry out the analysis described above, we first ex- 
amine our simulation trajectories to accumulate statistics 
for each time interval r, classifying type 2 particles ac- 
cording to the logarithm of their mean-square displace- 
ment during the observation period. This amounts to 
creating a histogram from the distributions shown in Fig- 
ures [5a] and I2bl assigning type 2 particles to "mobility 
bins" . Depending on the value of r, we use bin sizes for 
log 10 (<5r 2 ) in the range 0.1 - 0.2, which we find is suffi- 
ciently narrow to capture the shapes of the mobility dis- 
tributions, but coarse enough to allow for excellent statis- 
tical sampling. Following Ref. 46, we compute an average 
pair correlation function, g2j(r), between the type 2 par- 
ticles in a particular r-dependent mobility bin and all 
surrounding particles of type j. In determining 7)2j (r), 
the relevant pair separations are computed from config- 
urations sampled uniformly in time throughout the time 
interval of length r. In this work, we use the ~ overbar 
to denote a quantity that describes an average dynamic 
structure surrounding the type 2 particles belonging to a 
specific mobility bin for a given time interval r. 

In order to convert the T-dependent structural infor- 
mation contained in the g2j (r) into a number that char- 
acterizes the degree of average pair translational order 
that a type 2 particle (in given a mobility bin) experi- 
ences during the course of its displacement, we compute 
—82, which we define as 

P poo 

(1) 

Here, p — (iVj + N 2 )/V is the total number density, and 
Xj is the mole fraction of component j. This measure 
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is a dynamic generalization of a static structural metric, 

(2) 

— , which quantifies the contribution to the excess en- 
tropy of a mixture arising from equilibrium pair corre- 
lations involving particles of type 2. Our motivation for 
using — sjj in this study comes from (1) the earlier empir- 
ical observatio n 47 ' 48 that the long-time tracer diffusivity 
of species i in equilibrium mixtures scales in a simple way 

(2) 

with the static measure —s 2 and (2) the wider literature 
demonstrating that excess entropy captures many of the 
effects that temperature, density, and confinement have 
on the tr ansp ort coefficients of equilibrium fluids (see, 
e.g., Refs. [33M38T ) . Although we focus exclusively on the 
quantity — S2 in this paper to characterize dynamic struc- 
ture, we have found that other commonly used structural 
order metrics^ calculable from the r-averaged g2j (r) pro- 
duce qualitatively similar results. 
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FIG. 3: (Color online) Metric for the average structural order, 
—82, that type 2 particles of the binary mixture discussed in 
the text experience during an observation time r as a function 
of the logarithm of their mean-square displacement log 10 (<5r 2 ). 
Data is for packing fractions (a) ip — 0.57 and (b) ip = 0.582. 
The numbers in the figure correspond to the value of r for the 
curve of the same color. The size of the symbol is proportional 
to the fraction of particles in the "mobility bin" centered at 
that value of log 10 (Sr 2 ). In (b), the horizontal dashed line 
at Sr 2 = 0.56 represents the boundary between mobile and 
immobile particles (see Figure [2b] and text). 



To establish a baseline, we first examine the connec- 
tion between displacement and average structural order 
(during displacement) for type 2 particles in the system 
at (p = 0.57, a supercooled fluid state point for which pro- 
nounced dynamic heterogeneities have not yet emerged. 
Figure [3a] shows that, for all times spanning from the 
plateau region to the diffusive limit of the mean-square 
displacement curve in Figure [1] there is a clear negative 
correlation between the average structural order — sjj sur- 
rounding a particle during an observation window r and 
how far it moves in that time frame. As might be ex- 
pected, the short-time curves are considerably steeper 
than longer-time curves. That is, for shorter time inter- 
vals, larger displacements are, on average, accompanied 
by progressively more local structural disordering (i.e., 
weakening of the pair correlations associated with the 
tagged particle). At longer times approaching the dif- 
fusive limit, particles sample a broader distribution con- 
taining larger mean-square displacements, and the cor- 
relation between average structure and dynamics of the 
tagged particle is reduced. 
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FIG. 4: (Color online) (a) The value of s 2 for Sr 2 = 0.56, 
S2,c = S2(<5r 2 = 0.56), as a function r. (b) Reduced measure 
of structural order S2/s2,c(t) versus log 10 (5r 2 ) for ip = 0.582. 
The data and symbols are the same as in Figure [3b] 

The behavior is different, however, for the fluid at 
ip = 0.582 where the bifurcated displacement distribu- 
tion characteristic of strong dynamic heterogeneities is 
observed (see Figure [3b| . In particular, there is now a 
sharp change in the slope of the correlation between aver- 
age structure and dynamics when one compares "immo- 
bile" versus "mobile" particles. For immobile caged par- 
ticles, there is again a strong negative correlation between 
structural order and particle displacement; i.e., larger vi- 
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brational displacements are accompanied by increasingly 
weaker average pair correlations with neighboring parti- 
cles. Moreover, the immobile particles sample a struc- 
tural environment that is, on average, more ordered than 
mobile particles on the same time scale. However, within 
the class of mobile particles at a given time, the corre- 
lation between structural order and mobility is consider- 
ably weaker. 

In fact, if we focus on the structure and dynamics of 
mobile particles, a minor secondary effect in Figure I3bl 
is also apparent. Specifically, those particles with in- 
termediate displacements (of the order of a single parti- 
cle diameter) on a time scale r can have slightly more 
average structural order than those with either smaller 
or larger displacements. This feature is likely due to 
recaging events, — where a new coordination shell is tem- 
porarily formed around a particle that has traveled just 
far enough to "break free" from its original set of nearest 
neighbors. 

Another prominent feature of Figure [3b] is that the 
curves for different r have the same generic shape and 
appear merely shifted in terms of their average struc- 
tural order. In fact, Fig @}a) shows that the average 
structural order for particles with displacements at the 
boundary between immobile and mobile regions, i.e., 
~ S2,c = —S2{Sr 2 — 0.56), increases logarithmically with 
r under these conditions. This result makes intuitive 
sense — more structural order is, on average, expected 
to surround particles that take a longer period of time 
to exhibit the same value of mean-square displacement. 
Interestingly, as we show in Fig@|b), the r dependence 
of Figure I3bl is approximately removed altogether if one 
simply plots the average structural order of particles of a 
given mobility class normalized by its value at the bound- 
ary between mobile and immobile particles, S2/s2,C> ver- 
sus \og 10 (Sr 2 ) . In other words, the overall "scale" of the 
structural order surrounding particles has a simple r de- 
pendence for the intermediate times where dynamic het- 
erogeneities are present (500 < r < 8000). The relative 
differences between the structures surrounding particles 
in different mobility classes, on the other hand, show only 
a very weak dependence on r. 

Although we do not yet have a theoretical model for 
predicting the average pair correlations of particles with 
different mobilities, the logarithmic time dependencies il- 
lustrated Fig. 2] are perhaps not too surprising when one 
considers the following. Mean-square displacements of a 
tagged particle over a time r can be expressed 6D(t)t. 
Little is known about how D(t) relates to structure, 
but long-time self-diffusivities D(oo) exhibit an approx- 
imately exponential dependence on the two-body excess 
entrop y 34 i 47 ' 48 for dense fluids. Thus, to first approxima- 
tion, one might expect that particles which undergo the 
same mean-square displacement for different time scales 
t share a constant value of the product rexp[s2]. Al- 
though this crude argument is consistent with the data 
of Fig. 2] and has some intuitive appeal, it is far from rig- 
orous. More theoretical work will be needed to provide a 
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FIG. 5: (Color online) Average number of nearest neighbors, 
ntot, for type 2 particles of the binary mixture discussed in 
the text during an observation time r as a function of the 
logarithm of their mean-square displacement log 10 (<5r 2 ). Data 
is for packing fractions (a) <p — 0.57 and (b) <p = 0.582. The 
symbols are the same as those in Figure [3] The size of the 
symbols is proportional to the fraction of particles belonging 
to each mobility bin over time interval r. 



comprehensive understanding of the connections between 
dynamic structure and single-particle mobility observed 
here. 

Given the above results, another question one might 
logically ask is whether a simpler measure like co- 
ordination number, which roughly characterizes lo- 
cal density surrounding a particle, might provide the 
same qualitative information as —~S2- To test this 
idea, we have also collected statistics on n tot = 
47r J2j x jP Jo mm 3 r2 92j(r)dr, the average number of near- 
est neighbors surrounding type 2 particles (in a given 
mobility bin on time interval r). Here, ^minj is the lo- 
cation of the first minimum in g%j (r) . Figures [5a] and 
[b] display n tot as a function of \og 1Q (5r 2 ) for ip = 0.57 
and 0.582, respectively. In contrast to — s%, the average 
coordination number during the displacement does not 
provide a clear indication of the different structural envi- 
ronments that surround mobile versus immobile particles 
on intermediate time scales. 

To summarize, we have presented computer simula- 
tion data on a model supercooled fluid showing a corre- 
lation between average structure and mean-square dis- 
placements that is qualitatively different for immobile 
and mobile particles over times relevant for dynamic het- 
erogeneities. Interestingly, although the overall scale of 
the structural order depends logarithmically on observa- 
tion time, the relative differences in structural order be- 
tween mobile and immobile particles are largely indepen- 
dent of time. We provide a simple rationale for the afore- 
mentioned logarithmic time dependence, but more theo- 
retical work is needed to fully understand these trends. 
Finally, we show that the aforementioned structural dif- 
ferences between particles with different mobilities are 
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not reflected by the coordination number. 

In addition to further work on this problem in the 
arena of molecular simulation, it might be interesting to 
eventually explore whether the dynamical rules of kincti- 
cally constrained lattice models for supercooled liquids 40 
produce behavior qualitatively consistent with that of 
Fig-Elfrom our molecular dynamics simulations. It might 
also be interesting to compare these results to those 
from confocal microscopy experiments of dense colloidal 
suspensions ^ 1 ^ 39 
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